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罗 水 峡谷 -丘陵 复合 地 瑶 的 地 震 差 动 相干 困 数 
便 型 及 规律 特征 与 影 呈 
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摘 要 :构建 了 鹤 水 峡谷 - 诺 陵 复合 地 和 貌 的 地 震 差 动 模拟 所 需 的 相干 函数 模型 。 首 先 ,提出 并 推导 
得 到 了 考虑 土 体 滤波 效应 和 竖 向 距离 的 惟 水 峡谷 - 详 陵 复合 地 貌 相 干 函 数 模 型 表达 式 ; 探 讨 并 阐述 
了 相干 性 变化 的 规律 特征 及 其 物理 本 质 。 然 后 ,基于 新 相干 模型 和 传递 函数 构建 地 下 功率 谱 竹 阵 ， 

生成 地 下 多 点 地 震动 ;通过 功率 谱 及 相干 函数 拟 合 验 证 了 模拟 结果 的 科学 性 与 合理 性 ;对 比分 析 了 
本 研究 提出 的 相干 函数 模型 对 多 ee 变异 性 的 影响 。 最 后 ,以 一 座 连 续 刚 构 桥 为 合 ,探讨 
了 新 相干 函 数 模型 对 结构 地 震 差 动 响应 的 影响 。 结 果 表 明 : 基 于 新 相干 函数 模型 的 多 点 地 震动 相 
干 性 降低 ,各 点 地 震动 的 空间 a 加 剧 了 结构 各 支承 点 运动 的 不 同步 性 ,最 终 导 致 结 
构 地 震 差 动 响 应 增 大 。 
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A coherence function model for simulating multi-support seismic 
motions of canyon-hill composite topography and characteristics 


CHEN Xinyu ,LIU Guohuan ” ,FEI Qixiang 
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300350 Tianjin ,China;2. School of Civil Engineering, Tianjin University,300350 Tianjin, China) 


Abstract: This paper innovatively provides a coherence function model, specially adapted for water-covered 
canyon-hill composite topography and for simulating the multi-support seismic motions. First, a new 
coherence function model considering the soil filtering effect and vertical distance is proposed and derived. 
The characteristics of coherence and physical meaning of the proposed model are further investigated and 
summarized. Then, the multi-support underground seismic motions are generated by decomposing the 
underground power spectral density matrix , which ls formed by combining the coherence function based on 
the derived model with the transfer function. Subsequently, the reasonability and rationality of the 
simulated results are verified by comparing the simulated power spectrum and the coherence function to the 


theoretical ones. Finally ,the influence of the presented coherence function model on the differential effect 
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of multi-support seismic motions and the dynamic response of a continuous rigid frame bridge is specially 


analyzed , respectively. The results show that the coherence among the simulated variable seismic motions is 


relatively reduced ,but the variability of seismic motions at each foundation is increased , which leads to the 


enlargement of the structural responses. 


Key words: canyon-hill composite topography; coherence function; multi-support seismic motion; 


selsmic response 


大 量 研究 表明 场地 的 地 形 和 地 质 条 件 对 地 震动 
影响 显著 ,因而 山区 地 震动 与 平原 地 震动 在 特征 上 
表现 出 较 大 差异 。 周 田 字 等 … 利 用 当前 较 具 代表 性 
的 地 形 效 应 量化 表征 模型 ,初步 分 析 了 川 汗 山 区 高 
频 地 震动 的 地 形 效应 ; 崔 殖 等 ”从 幅 值 持 时 和 频谱 
3 个 角度 探讨 了 高 山峡 谷地 区 洞 室 群 场 址 的 地 震动 
特征 ; 王 运 生 等 中 基于 实测 资料 对 丘陵 地 貌 区 斜坡 
地 震动 响应 特征 开展 了 研究 ,证 实 了 斜坡 高 位 放大 
效应 , 且 表 明 坡 顶 的 高 频 滤波 作用 更 为 明显 。 受 地 
形 效应 , 行 波 效应 和 部 分 相干 效应 的 影响 ,地 震动 表 
现 出 明显 的 空间 变异 性 , 由 于 缺乏 各 类 地 形 条 件 下 
的 多 点 地 震动 实测 资料 , 故 研究 考虑 地 形 效应 的 多 
点 地 震动 人 工 模拟 方法 对 山区 结构 尤其 是 大 坝 .大 
跨 桥梁 和 长 距离 输电 塔 等 长 大 结构 的 抗震 分 析 具 有 
重要 意义 。 通 常 多 点 地 震动 模拟 需要 通过 目标 功率 
谱 , 传 递 函 数 和 相干 函数 构建 功率 谱 和 矩阵 ,基于 此 研 
究 框架 ,WU 等 中 利用 水 平地 表 到 峡谷 表面 的 传递 
函数 模拟 了 V 形 峡 谷 表面 的 多 点 地 震动 ,其 相干 函 
数 采 用 了 目标 点 的 水 平 投影 点 来 近似 计算 ; LIU 
等 提出 了 峡谷 相干 函数 并 进一步 计算 得 到 地 下 
相干 函数 ,利用 水 平地 表 到 峡谷 表面 再 到 地 下 的 两 
步 传递 函数 法 模拟 了 形 峡 谷 场地 的 地 下 多 点 地 震 
动 ; 柳 国 环 等 ”分 别 获得 了 分 层 圆 弧 峡 谷 在 SH 波 
和 SV 波 入 射 下 的 多 点 地 震动 ,进而 获得 了 V 形 峡 
谷 的 地 下 多 点 地 震动 ,并 分 析 了 其 上 桥梁 在 超大 震 
下 的 破坏 模式 ; 何 颖 等 提出 了 一 种 考虑 土 体 非 线 
性 特性 的 沉积 河谷 多 点 地 震动 模拟 方法 ,并 以 V 形 
半 满 河谷 为 例 验 证 了 其 合理 性 。 

以 上 有 关 多 点 地 震动 模拟 的 研究 , 均 针对 单一 
峡谷 地 形 ,而 在 山区 往往 存在 相距 不 远 的 多 个 峡谷 ， 
各 峡谷 间 可 视 作 丘陵 ,形成 一 种 峡谷 -丘陵 复合 地 
貌 。 学 者 们 针对 该 类 复合 地 貌 开展 了 一 些 研究 : 李 
郑 梁 等 “采用 边界 积分 方程 法 得 到 了 浅 切 割 高 山 
峡谷 复合 地 貌 的 位 移 响应 ,并 探讨 了 入 射 波 类 型 及 
频率 、 入 射 角 、 山 体 和 峡谷 的 几何 参数 对 该 复合 地 貌 


地 震动 放大 效应 的 影响 ; 张 玉 人 敏 等 2 拟 合 了 高 山峡 
谷地 区 地 下 洞 室 群 非 平稳 人 工地 震动 ,但 其 模拟 的 
并 非 空间 相关 的 多 点 地 震动 。 故 尚 需 发 展 一 种 峡 
谷 -丘陵 复合 地 貌 地 下 多 点 地 震动 的 有 效 模拟 方法 ， 
其 关键 在 于 获得 适用 于 该 复合 地 貌 的 相干 函数 。 在 
V 形 峡 谷 多 点 地 震动 模拟 中 ,峡谷 表面 两 点 的 相干 
函数 采用 了 其 水 平地 表 投 影 点 的 相干 函数 近似 痊 
代 , 而 对 于 峡谷 - 斤 陵 复 合 地 貌 , 需 考 虑 括 陵 的 滤波 
作用 ,因而 现 有 的 相干 函数 计算 方法 不 再 适用 。LIU 
等 ”提出 了 一 种 通过 位 移 函 数 的 传 里 叶 谱 计算 峡谷 
相干 函数 的 方法 ,但 其 依赖 于 场地 的 散射 理论 解 , 且 
无 实测 资料 提供 依据 。 因 此 尚 需 发 展 一 种 基于 经 验 
相干 函数 模型 ,并 需 考 虑 峡谷 -丘陵 复合 地 貌 特征 的 
新 相干 函数 模型 。 

本 研究 以 覆 水 双 V 形 峡谷 为 例 ,首先 推导 了 和 覆 
水 峡谷 -丘陵 复合 地 貌 相干 函数 的 显 式 表达 式 , 对 该 
相干 函数 模型 的 优势 加 以 分 析 并 探讨 了 其 相干 性 的 
变化 规律 和 物理 本 质 ;进而 模拟 得 到 了 基于 该 相干 
函数 模型 的 地 下 多 点 地 震动 ,通过 功率 谱 和 相干 函 
数 拟 合 验证 了 所 模拟 地 震动 的 谱 兼 容 性 和 空间 相关 
性 ;对 比分 析 了 新 相干 函数 模型 对 多 点 地 震动 空间 
变异 性 的 影响 ;最 后 ,以 一 座 连续 刚 构 桥 为 例 , 探 讨 
了 新 相干 函数 对 结构 差 动 响应 的 影响 。 本 研究 可 为 
我 国 西南 部 峡谷 -丘陵 地 貌 区 的 多 点 地 震动 模拟 和 
相应 长 大 结构 的 抗震 分 析 提供 参考 。 


1 覆 水 峡谷 -丘陵 复合 地 貌 相干 函数 


1.1 推导 过 程 


场地 模型 如 图 1 所 示 , 要 得 到 点 i 和 n 的 相干 函 
数 , 取 位 于 峡谷 另 一 侧 且 与 点 i 高度 相同 的 点 j, 则 
点 7 的 振动 可 以 用 具有 相位 差 的 点 i 的 振动 来 表 
示 四 (注意 到 7， = 7 ,hh = hs ) , 即 
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(w) AW, (w) p 
7 二 Ary) i (++ Li 
油 Ti 天” wi 7 


A | 

ey (1 + Ar7,) (1) 
其 中 :4 代表 谐 波 运动 的 振幅 ; 7; 和 7 为 常数 ; 7 = 
k*z/w = z* Vp/C” ,p 代表 目标 点 处 土壤 密度 ,z 为 
目标 点 路 地 表 的 深度 ,CG” 和 "分别 为 复 波 数 和 复 
剪 切 模 量 ; p, 为 水 的 密度 ; 有 为 水 的 复 体积 模 量 ; 
h,, 和 刀 , 分 别 为 i 、j 两 点 距 履 水 水 面 的 深度 。 


和 


图 1 场地 示意 图 


Fig.1 Schematic diagram of site 
根据 随机 振动 理论 , i 、j 两 点 的 互相 关 亲 数 可 
以 表示 为 
民 (Ari) = ELi®® (六 (GTAri)] 


》 
se AR 


宣 儿 R,( Ar,) (2) 
其 中 R;(A7,) 为 点 的 目 相 关 困 数 。 
进而 i 、7 两 点 的 互 功 率 详 密 度 函 数 可 表示 为 
Si(o) = | Ry(Ari)e "id(Ar,) 


A -ik 
局 本 Ri(AT,)e "dd(Ar,) 
(w) 
Sa(®) (3) 
4 
由 相干 函数 的 定义 , S;(w) 还 可 以 表示 为 
Si(o) = yo) VSi(w)S;(w) (4) 


联 立 式 (3) 和 (4) ,得 到 
A 
V Sj;(w) = 汉 AAS) (5) 
A ysl(w) 
相似 地 , m 、n 两 点 间 也 存在 如 下 关系 , 即 


A yi ©) 
ld VS 0) (6) 


(Cn) 
4 


假设 点 7 和 m 间 的 土 体 为 单 自由 度 结构 ,地 震 
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波 由 点 7 问 m 传播 ,点 7 处 的 地 震动 加 速度 aj(t) 经 
过 土 体 滤波 后 ,得 到 点 m 处 的 相对 位 移 、 速 度 和 加 速 
度 时 程 分 别 为 m(?) 、m'(t) 和 到 , 则 土 体 的 振 
动 方程 为 

m(t) +2om (t) + wm(t) = 一 a(t) (7) 
式 中 4 和 wo 分 别 为 土 体 的 阻尼 比 和 目 振 频率 。 

对 式 (7) 两 侧 进行 傅 里 叶 变 换 , 得 到 相对 位 移 、 
速度 和 加 速度 的 傅 里 叶 详 分 别 为 
— 1 


m(w) = er (8) 
人 = 了 5 
m"(w) = 和 ia(w) (10) 


2 
一 ww +1 .25auuw + wo 


点 m 处 的 绝对 加 速度 为 


an = mt) + a(t) =- 200m (1) 一 wom(t) 
(11) 

对 陈 (11) 两 侧 进行 传 里 叶 变 换 得 到 
a (0) =m (oO) ta(w) = -200m'(w) -wm(w) 
(12) 


将 式 (8) 和 (9) 代 入 式 (12) ,得 到 
i + 2 + wo 
二 < 27W0w + oe 
= H;,(iw)a,(w) (13) 
进而 ,点 j 和 m 的 自 功率 谱 间 有 如 下 关系 “ 
Sia(®) = Hi; (iw)H; (iw)S;(w) (14) 
将 式 (5) 和 (6) 代 入 式 (14) 得 到 


an,(wW) = 


1 

VH, (iw) H(iw) 

Sj(w) 

ON Sinl@) 

A 1 Si(w) 
Ad 


进一步 整理 得 到 
A SU) H(iw) 
Cw) i 一 Im (@w 0 (16 
最 终 i 、n 两 点 的 互 功率 谱 密 度 孔 数 可 以 表示 
为 (注音 到 此 时 7; 关 7, ,hh 2 有, ) 
Lis Dw Dw 
Si(iw) = Te VV Sw) (17) 
(w) 


由 式 (16) 和 (17), 可 得 i、n 两 点 相干 也 数 为 
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yi iw) yun(iw) =exp[i0;,,(w)| x 
so Ri 
VSi(w)S,,(w) 其 中 
A (n/a Si(w) en ho | 
AD® wr Bm) So 0,,(w) = arctan RN (21 ) 
H. io(Tn+ Mpw/ Kv hwn-Ti- Mpw/ Kv hwi) 意 日 
= ye) ee 
hiw) 征 目 标点 的 竖 向 高 度 差 对 相干 函数 的 影响 。 说 明 本 


(18) 

其 中 yj(w) 和 Yi,(w) 采用 Hao 模型 计算 。 

式 (18) 即 本 研究 提出 的 适用 于 上 覆 水 峡谷 -丘陵 
复合 地 貌 的 地 表 相 干 遇 数 计算 方法 。 

在 以 往 研 究 中 ,考虑 到 了 波 可 在 水 中 传播 ,从 随 
机 振动 理论 和 一 维 波动 理论 出 发 ,推导 了 和 窗 水 效应 
的 地 下 相干 函数 ,其 具体 表达 式 为 

y,, (iw)e™ Ge 


(目标 点 位 于 无 覆 水 区 ) 


p 
(C77y) (TT + x (hnn’ —hiz') ] 
W 


yi (IW) = 
Yn(iw)e™ 
(目标 点 位 于 覆 水 区 ) 
(19) 
由 式 (18) 计算 得 到 地 表 点 间 的 相干 函数 
ya(iw) 后 ,由 式 (19) 即 得 到 相应 地 下 点 间 的 相干 函 
数 yim (10w) O 
1.2 改进 前 后 相干 函数 的 对 比 
式 (18) 在 峡谷 部 位 相干 函数 的 基础 上 , 乘 以 因子 
VH "(iw) /Bn iw) 和 和 
到 峡谷 - 顾 陵 复 合 地 貌 相 干 函 数 , 其 中 因子 
VD,(iw)/H;(iw) 的 物理 意义 是 表征 丘陵 土 体 的 滤 
波 效 应 ,该 因子 的 模 为 1 ,表明 丘陵 部 位 的 土 体 仅 影 响 


相位 ,不 对 捞 清 相干 冰 数 产生 影响 ,该 因子 可 以 写 为 
相位 角 的 形式 。 


研究 提出 的 相干 函数 ,不 仅 考 虑 了 目标 点 间 的 水 平 距 
离 , 还 计 入 了 竖 向 距离 的 影响 。 
为 探究 改进 后 相干 性 的 变化 规律 ,选取 如 图 2 所 
示 的 目标 点 ,各 点 坐标 已 在 图 中 注 明 。 图 3 给 出 了 改 
进 前 后 地 表 点 4D .4C 和 BD 间 的 相干 函数 ,其 中 改进 


前 的 相干 也 数 即 由 目标 点 的 水 平 投影 点 计算 ,改进 后 
的 相干 函数 按 式 (18 ) 计算。 
区 人 Eo 90,80) G080) a 80) 


> 人 
! | B(—90,0) 0 1 从 
业 上 由 由 
A!’ 万 你 2 


A | 
(~120,40) 


下 二 


图 2 日 标 点 示意 轿 
Fig.2 Schematic diagram of target points 
由 图 3 可 知 ,与 水 平 投影 点 计算 得 到 的 相干 也 数 
相 比 ,本 人 研究 提出 的 覆 水 峡谷 -于 陵 复合 地 貌 相干 也 
数 整体 上 表现 出 相干 性 减弱 ,但 在 频率 很 小 时 ,本 人 研 
究 相 干 函 数 的 相干 性 略 有 增强 。 其 物理 本 质 在 于 :在 
相干 函数 中 引入 竖 癌 距离 的 影响 后 ,各 点 地 震波 传播 
路 径 的 差异 增 大 ,进而 传播 过 程 中 地 震波 的 幅 值 与 相 


位 角 产 生变 异 的 概率 随 之 增 大 。 同 时 频率 越 高 ,传播 
过 程 中 振动 次 数 越 多 , 幅 值 与 相位 角 产 生变 异 的 概率 
也 越 大 。 当 频率 很 小 时 , 竖 癌 距离 对 相干 性 的 影响 可 
以 忽略 ,此 时 相干 性 反而 略 有 增强 。 


0 2.5 5.0 5 00 2s 0 D3 0 75 7100 1125 0 2.5 50 75 100 195 
频率 /Hz 频率 /Hz 频率 /Hz 
(a) .47 地 表 相 干 函数 (b) 4C 地 表 相 于 函数 (c) 3D 地 表 相 干 函数 


图 3 改进 前 后 相干 性 的 变化 


Fig.3 ”Comparison of coherence before and after Improvement 
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2 和 窄 水 峡谷 -丘陵 复合 地 貌 多 扣 地 震动 模拟 


2.1 基于 本 研究 模型 的 多 点 地 震动 模拟 


依据 地 震动 模拟 理论 ,模拟 了 和 宪 水 峡谷 - 左 院 复 
合 地 貌 的 多 点 地 震动 ,目标 点 的 位 置 及 坐标 见 图 2。 
地 震动 持 时 取 为 20. 48 s, 选 取 Clough-Penzien 功率 
谱 模 型 “作为 水 平地 表 的 目标 功率 谱 , 其 参数 取 为 
So = 0.0069,w, =57 rad/s, w=0. lw,,é, =é,= 
0.6。Hao 相干 模型 ”参数 采用 SMART-1 台 阵 45 
条 地 震 记 录 的 统计 参数 “ ,其 值 B,、a.b 和 * 分 别 取 
为 1.109 x10”“、3.853 x10、 -1.811 x10- 和 
1. 177 x10“。 选 取 Jennings 包 络 函数 "对 生成 的 
地 震动 进行 调幅 ,其 参数 4 以 及 c 分 别 取 为 1.2.9 


2 地 表 4 点 地 震动 
一 一 地 下 4' 反 地 震动 


PGA,,.=1.008 
Sy On 二 0 4 27 8 


sl 
9 
= 0 
芍 
| 
a PGA =0.863 
OwF—0.209 
2 
. 10 [3 20 20.48 
时 间 /s 
必 POA 1333 地 表 C 点 地 震动 
Ow#—0.323 只 一 一 地 FC' 点 地 震动 
号 
~ 
襄 
芍 
be 
上 
10 ]5 20 20 48 
时 间 /s 


和 0.155。V 形 峡 谷 半 和 覆 水 ,水 的 体积 模 量 取 2.2 x 
10”Pa, 密 度 取 1 000 kg/m , 土 体 物理 参数 见 表 1。 
图 4 给 出 了 各 目标 点 的 模拟 加 速度 时 程 ,图 中 标注 
的 PCA 代表 加 速度 峰值 ,o 代表 所 模拟 加 速度 时 程 
的 标准 差 。 结 果 表 明 , 各 地 下 目标 点 的 标准 差 均 小 
于 相应 的 地 表 目 标点 ,这 是 由 土 体 的 滤波 作用 和 放 
大 效应 导致 的 。 此 外 ,峡谷 地 形 表现 出 了 明显 的 放 
大 效应 ,接近 谷 顶 的 C.D 两 点 地 震动 大 于 接近 谷底 
的 4.B 两 点 。 


表 1 各 层 土 体 物理 参数 


Tab.1 Physical parameters of soil 
层 ”密度 p/(kg."m ) 剪 切 模 量 G/MPa 阻尼 比 5 
5 2 500 380 0.05 
d 2 650 320 0.05 
2r =----- 地 表 8B 点 地 震动 


一 一 地 FB' 点 地 震动 


-| 
目 
= 
芍 
| PGAww-0.965 
“APGA, .=1.179 Ow#™0.2069 
Our=0.240 
2 
, 5 10 15 20 20.48 
时 | 间 /s 
7 地 表 几 点 地 震动 
一 一 地 下 DD' 点 地 震动 
| 
8 
二 0 
EX | 
网 ' | 
Ry , FGA 1.307 
PGA，=1.120 | 一 1 re 
ou -0.253 
,|L_ 0-0.299 
y 5 10 15 20 20.48 


时 间 /s 


图 4 地 表 及 地 下 目标 点 模拟 加 速度 时 程 


Fig.4 Ground and underground acceleration time histories of target points 


2.2 功率 谱 及 相干 函数 拟 合 
图 5 给 出 了 各 点 加 速度 时 程 的 频谱 与 目标 功率 


0.04 


a 地 表 模 拟 谱 
区 一 地表 目标 谱 
名 0.03 

0.02 

9 

= 

[a 


频率 /Hz 
(a) 地 表 有 点 功率 谱 


投稿 网 站 :http://cjam. xjtu. edu. cn 


谱 的 对 比 ,限于 篇 幅 , 仅 以 点 B(B') 和 C(C') 为 例 。 
结果 表明 ,所 模拟 加 速度 时 程 的 功率 谱 与 目标 功率 


\ 少 世 | 人 和 
谱 拟 合 良 好 。 
oD 
a 地 表 模拟 谱 
本 一 “地表 目标 谱 
号 0.03 
0.02 
5 
vl 
户 3 


频率 /Hz 
(b) 地 表 C 点 功率 谱 
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第 
0.04 hn 0.04 搬 上 二 | 让 
Ws 地 表 模 拟 详 | 2 i 
' 一 一 地表 目 标 并 . 
0 - 二 0.03 
2 0.02 


PSD/((m 
Det 
已 


0 i 
频率 /Hz 
(d) 地 下 C ' 点 功率 谱 


0 2.5 0 a0 1 
频率 /Hz 
(c) 地 下 B' 点 功率 谱 


图 5 目标 点 功率 谱 拟 合 


Fig.3 PSD fitting of simulated results with target PSD 
图 6 给 出 了 相干 函数 的 拟 合 情 况 ,限于 篇 幅 , 以 。 但 总 体 上 模拟 值 仍 围绕 目标 值 变 化 , 故 所 模拟 地 震 
生 , 证 明了 本 人 研究 所 模 


AB4'B'、BD 和 B'D' 间 的 相干 孙 数 为 例 。 绪 果 表 ” 动 表现 出 了 良好 的 空间 相关 + 
明 , 尽 管 相干 函数 的 模拟 值 出 现 了 一 定 程 度 的 波动 ， ” 拟 地 震动 的 科学 性 与 合理 性 


- -~ -~ AB 模拟 相干 函数 一 ……… A4'B' 模 拟 相 干冰 数 -= -BD 模拟 相干 函数 ……… B'D' 模 拟 相 于 通 数 
Le 4B 目标 相干 函数 一 一 4'B' 目 标 相干 函数 1<[ ----- BD 目标 相干 函数 一 一 B'D' 目 标 相干 函数 


5.0 3 10.0 此 5 


0 2 5.0 VS 10.0 |25 0 2 
频率 /Hz 


频率 /Hz 
图 6 ”地表 及 地 下 相干 函数 拟 合 


Fig.6 Fitting of ground and underground coherence functions 
ee 示 , 左 右 两 侧 分 别 给 出 了 改进 前 后 点 4' 和 D' .4' 和 
i 5 以 及 忆 和 六 地 筑 动 人 的 对 比 ,其 中 改进 的 
地 震动 是 采用 目标 点 水 平 投影 点 的 相干 函数 模拟 得 
到 的 ;改进 后 的 地 震动 是 采用 本 人 研究 提出 的 峡谷 - 丘 


3 间 变 异性 的 变化 规律 。 如 图 7 所 陵 复 合 地 貌 相干 函数 模拟 得 到 的 。 
27 | 
一 一 A' 点 地 震动 


各 点 地 震动 的 空间 变异 性 引起 差 动 效应 是 多 点 
震 输 入 与 一 致 输入 的 核心 区 别 , 夏 有 必要 人 研究 改 
进 后 多 点 地 震动 空 


2 D' 点 地 震动 Ge 
2 i 一 一 4’ 点 地 震动 ~ IO0%211% 
: 


| 5 
E E 
Sg Sa 
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MPGA， -0997 wy 0 6 POA, =0.863 -pt hp 1307 
0, =0.215 0,, =0.209 0p =0.253 
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0 5 10 15 20 0 5 10 15 20 
时 间 /s 时 间 /s 
(a) 改进 前 4' 和 D’ 间 (b) 改进 后 4' 和 D’ 间 
2 C' 点 地 震动 2 C' 点 地 震 亏 
一 才 / 二 攻 屠 元 Oo — O04 一 一 震 5 
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BO 
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国医 
0 5 10 15 20 0 5 10 15 20 


时 间 /s 
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(c) 改进 前 4 和 C ' 间 


投稿 网 站 :http://cjam. xjtu. edu. cn ” 微 信 公众 号 :应 用 力学 学 报 


398 应 用 力学 
WE J DD' 点 地 震动 
PGA, =1.119 一 一 B' 点 地 震动 
了 0,, =0.238 
| 
人 ' 
虽 A 
0 
be | b 
雹 ! 
| 0 0 
PGA,, =0.859 “X100%= -0.8% 
on =0.236 4 
“0 5 10 15 20 
时 间 /s 


(e) 改进 前 3' 和 D' 间 


.Ee DD' 点 地 震动 
2 = fh 元 
GX100%=5.4% 一 一 B' 点 地 震动 
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SEE0 
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0,, =0.240 
20 
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10 
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(改进 后 B' 和 DD' 间 


图 7 改进 前 后 各 点 地 震动 差 值 的 对 比 


Fig.7 Comparison of seismic motion variability based on before-and-after Improvement model 


由 图 7 可 知 , 本 研究 提出 的 相干 函数 增 大 了 各 
点 地 震动 的 幅 值 差异 ,进而 加 剧 了 多 点 地 震动 的 空 
间 变 异性。 改进 前 点 4' 和 D’'.A4' 和 C' 以 及 B' 和 D' 
的 地 震动 标准 差分 别 相 差 9.8% 、17.7% 和 0. 8%， 
改进 后 这 一 数值 变 为 21. 1% 、25. 8% 和 5.4%。 其 
原因 在 于 ,地 震动 可 以 看 作 一 系列 谐 波 运动 的 鳃 加 ， 
本 人 研 究 提出 的 相干 函数 考虑 竖 癌 距离 , 谐 流 传播 路 
径 差 异 增 大 ,相干 性 降低 ,两 点 同 频 谐 波 分 量 的 幅 值 
和 相位 角 差 异 增 大 ,最 终 导 致 两 点 地 震动 幅 值 和 标 
准 差 的 差异 增 大 。 


3 ” 算 例 分 析 


本 章 以 一 座 连 续 刚 构 桥 为 例 , 探 守 峡谷 -丘陵 复 
合 地 貌 相 干 函 数 对 多 点 地 震动 进而 对 结构 差 动 啊 应 
的 有 影响。 为 了 独立 考察 单一 因 系 (地 震动 ) 不同 输 


三 
页 0 


1 号 跨 


三 
页 和 


入 对 结构 反应 的 影响 ,从 而 充分 体现 以 上 研究 的 必 
要 性 和 重要 性 ,由 此 必须 避免 材料 进入 塑性 这 一 影 
响 因素 对 结构 反应 产生 影响 的 混淆 。 因 此 ,本 研究 


对 结构 只 进行 弹性 时 程 分 析 。 


= 
页 和 


3.1 工程 概况 与 有 限 元 模型 


如 图 8 所 示 , 某 5 监 钢筋 混凝土 连续 刚 构 桥 路 
越 峡 谷 - 撕 陵 复 合 地 貌 , 半 禾 水 V 形 峡 谷 宽 120 m , 深 
80 m, 土 体 物理 参数 在 图 8 中 注 明 。 刚 构 桥 全 长 
300 m, 跨 径 组 合 为 35 m +75 m+80 m+75 m+35m, 
采用 变 截 面 箱 染 , 箱 梁 顶 部 设 2% 的 单身 横 坡 , 墩 顶 
及 跨 中 截面 尺寸 见 图 9。1 号 和 4 号 墩 高 38 m, 墩 底 
埋 深 21 m,2 号 和 3 号 墩 高 19 m, 壤 底 埋 深 17 m, 桥 
墩 采 用 8.5 m x4 m 算 形 断面 。 混 凝 土 采 用 C50, 纵 
筋 和 链 筋 均 采 用 HRB400 。 


密度 p=2 500 kg/m 
剪 切 模 量 G=380 MPa 
泊 松 比 y=0.22 


密度 Po-2 650 kg/m 
甬 切 模 量 Ce-320 MPa 
泊 松 比 盖 0.22 


单位 : m 


图 8 某 连 续 刚 构 桥 立 面 图 


Fig.8 Elevation of a continuous rigid frame bridge 
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1 单位 : em 


图 9 主 架 墩 顶 ( 左 ) 及 路 中 ( 右 ) 截 面 
Fig.9 Section of beam at pier top (left) and midspan (right ) 
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采用 ABAQUS 有 限 元 软件 建立 了 桥梁 的 数值 
计算 模型 ,如 图 10 和 11 所 示 。 混 凝 土 结构 和 土 体 
均 采 用 C3D8I 实体 单元 , 共 21 692 个 ;钢筋 采用 
T3D2 机 架 单 元 , 共 7792 个 。 墩 梁 绑 定 约束 ,钢筋 般 
人 混凝土 ,由 于 本 研究 不 关注 桩 土 相互 作用 , 桩 土 采 
用 共 节 点 简化 处 理 。 深 体 两 侧 只 施加 竖 向 位 移 约束 
及 绕 桥 回转 动 约束 (U3 = UR1 = UR2 =0) ,在 所 截取 
地 基 的 底 边 界 和 侧 边 界 施加 接地 的 弹 先 和 阻尼 器 以 
建立 黏 弹性 人 工 边 界 ,模拟 地 基 的 辐射 阻尼 。 黏 弹性 
人 工 边界 的 刚度 系数 和 阻尼 系数 按 文献 [16] 计 算 。 


地基 放大 图 


图 10 桥梁 有 限 元 模型 
Fig. 10 Finite element model of bridge 


1 提取 梁 中 轴线 处 下 部 纵 盘 
I 各 截面 最 大 正 应 力 ， 如 网 1S(b) 


图 11 钢筋 骨架 


Finite element model of reinforcement skeleton 


Fig. 11 
3.2 ”本 研究 提出 的 相干 函数 对 结构 响应 的 影响 


为 进一步 探究 本 研究 提出 的 履 水 峡谷 -丘陵 复 
合 地 貌 相干 函数 对 结构 响应 的 有 影响, 分别 采 用 改进 
前 的 相干 函数 和 本 研究 提出 的 相干 函数 模拟 了 各 桥 
墩 支 承 点 的 地 下 多 点 地 震动 ,如 图 12 所 示 , 计 算 并 
对 比 了 两 组 地 震动 激励 下 刚 构 桥 的 动力 响应 。 选 取 
2.3 和 4 号 跨 的 跨 中 和 两 端 ,以 及 1.3 两 桥墩 的 墩 顶 
和 墩 底 作 为 观测 点 。 输 入 前 首先 对 各 组 地 震动 进行 
归 一 化 处 理 , 即 对 每 一 组 地 震动 数据 除 以 该 组 数据 
中 的 最 大 值 ,再 对 归 一 化 后 的 数据 分 别 按 加 速度 峰 
值 为 0.2g 和 0.4g 进行 输入 ( 为 重力 加 速度 ) 。 所 
模拟 地 震动 入 射 角 为 30° ,将 其 分 解 为 顺 桥 向 和 竖 回 


提取 染 中 轴线 处 上 部 纵 篇 轴 癌 的 
各 截 5 


面 最 大 正 应 力 ， 如 图 15(a) 


有 限 元 模型 


从 人 。 参 照 文献 117 | 将 地 震动 转化 为 下 接 作 用 于 
边界 市 点 上 的 等 效 集 载 来 实现 波动 输入 。 

表 2 给 出 了 改进 前 后 刚 构 桥 各 观测 点 最 大 内 力 
的 对 比 。 表 3 给 出 了 改进 前 后 各 观测 点 最 大 位 移 的 
对 比 ,图 13 和 14 分 别 为 改进 前 后 0.2g 地 震动 峰值 
下 1 号 桥墩 底部 和 顶部 的 顺 桥 向 位 移 时 程 和 加 速度 
时 程 的 对 比 。 结 有 果 表 明 , 使 用 峡谷 - 顾 陵 复合 地 貌 相 
干 孔 数 后 ,本 研究 刚 构 桥 的 差 动 响应 普遍 增 大 ,一 方 
面 是 由 于 改进 后 2、3 号 桥墩 处 地 震动 显著 增 大 ,为 
一 方面 是 由 于 各 点 地 震动 的 空间 变异 性 增强 ,加 剧 
了 各 文 承 点 运动 的 不 同步 性 。 图 12 中 改进 前 后 的 
地 震动 最 大 仅 相 差 8.3% (2 号 桥墩 处 ) ,而 结构 响应 
的 最 大 增幅 达到 了 118.3% (0. 4g 峰值 激励 下 3 号 
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跨 跨 中 竖 向 位 移 ) ,由 此 表明 ,多 点 输入 时 结构 差 动 进 前 的 相干 国 数 将 低估 结构 啊 应 。 


多 ee 
PEGA 10190 改进 前 
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图 12 改进 前 后 输入 的 多 点 地 震动 
Fig. 12 Multi-support seismic excitations based on before-and-after improvement model 
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图 13 顺 桥 向 位 移 时 程 对 比 
Fig. 13 Comparison of along-bridge displacement time history 
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图 14 顺 桥 向 加 速度 时 程 对 比 
Fig.14 Comparison of along-bridge acceleration time history 
图 15 给 出 了 改进 前 后 0.2g 地 震动 峰值 下 沿 染 ”” 进 前 后 0.2g 地 震动 峰值 下 的 墩 底 剪 力 时 程 对 比 和 
内 纵 筋 轴 向 的 各 截面 最 大 正 应 力 对 比 , 所 提取 钢筋 ” 沿 桥墩 高 度 的 最 大 剪 力 对 比 , 限 于 篇 幅 , 仪 以 2 号 和 
的 位 置 在 图 11 中 注 明 。 图 16 和 17 分 别 给 出 了 改 ”3 号 桥墩 为 例 。 结 采 表 明 ,使 用 本 研究 提出 的 相干 
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阴 数 模型 后 ,2 号 跋 和 3 号 蜂 处 的 钢筋 截面 最 大 正 ”本 人 研究 提出 的 覆 水 峡谷 -丘陵 复合 地 貌 相 干 郴 数 对 
应 力 以 及 桥墩 的 最 大 剪 力 显著 增 大 ,进一步 证 明了 实际 工程 具有 不 可 忽视 的 意义 

表 2 改进 前 后 的 最 大 内 力 对 比 


Tab.2 “ Comparison between maximum internal force obtained from before-and-after Improvement model 


0.2g 0.4g 
剪 力 / 弯 和 矩 / 剪 力 / 弯 矩 / 
观测 点 Fi-F, M| -Mo 下 一 有 1 - M, 
(107kN) 本 (102kN . m) 去 /% (107kN) (10°kN . m) Dn [OO 
0 Fn Mo MI i Fn Mo MI 
2 跨 左 端 240.7 357.1 48.4 386.7 810.5 109.6 
2 跨 右 端 ”394.4 635.7 61.2 607.5 1091.9 79.7 
3 跨 左 端 332.7 543.9 63.5 521.0 944.3 81.2 
3 跨 右 端 ”335.6 363.8 8.4 526.8 656.1 24.5 
4 跨 左 端 241.3 341.9 41.7 621.4 652.8 5.1 
4 跨 右 端 501.2 276.3 -44.9 908.1 457.6 -49.6 
2 跨 跨 中 2 858.7 2 882.7 0.8 4965.0 4758.6 -4.2 
3 跨 跨 中 1513.6 2207.4 45.8 2933.8 4501.4 53.4 
4 跨 蜂 中 2349.8 1807.8 -23.1 4682.3 3739.1 -20.1 
1 墩 顶部 ”243.8 391.0 60.4 10535.6 15502.9 47.1 442.4 827.2 87.0 19181.0 32905.2 71.6 
1 墩 底部 258.5 399.0 54.4 1090.6 730.0 -33.1 471.7 843.7 78.9 2006.2 1630.4 -18.7 
3 墩 顶部 278.1 344.4 23.8 7208.7 8166.2 13.3 578.5 687.8 18.9 15022.7 15728.1 4.7 
3 墩 底部 ”292.2 352.6 20.7 1798.9 2427.8 35.0 606.7 693.3 14.3 3781.1 4668.0 23.5 
注 :代表 改进 前 的 剪 力 ;下 代表 改进 后 的 剪 力 ;Mo 代表 改进 前 的 弯 矩 ;Mi 代表 改进 后 的 弯 矩 。 
表 3 改进 前 后 的 最 大 位 移 对 比 
Tab.3 Comparison between maximum displacement obtained from before-and-after improvement model 
0.2g 0.4g 
顺 桥 向 位 移 /m 竖 向 位 移 ]m 顺 桥 向 位 移 /m 坚 向 位 移 ]m 
观测 点 六- 二 ZE Yi 一 芒 Zi -Zo 
yy yy WY 有 站 
2 跨 跨 中 0.048 0.087 81.3 0.122 0.155 27.0 0.095 0.174 83.2 0.221 “0.331 49.8 
3 跨 跨 中 0.047 0.092 95.7 0.096 0.183 90.6 0.095 0.185 94.7 0.175 “0.382 118.3 
4 跨 跨 中 0.048 0.091 89.6 0.081 0.118 45.7 0.096 0.182 89.6 0.141 ”0.258 83.0 
1 墩 顶部 0.057 0.100 75.4 0.087 ”0.090 3.4 0.112 0.202 80.4 0.160 0.195 21.9 
1 墩 底部 0.057 0.086 50.9 0.087 “0.091 4.6 0.112 ”0.174 55.4 0.160 0.196 22.5 
3 墩 顶 部 0.045 0.086 91.1 0.087 ”0.149 71.3 0.090 ”0.172 91.1 0.157 0.314 100.0 
3 墩 底部 0.042 0.085 102.4 0.087 ”0.149 71.3 0.084 0.170 102.4 0.157 0.314 100.0 


注 : 刀 代表 改进 前 顺 桥 向 位 移 ; 六 代表 改进 后 顺 桥 向 位 移 ;2Z0 代 表 改 进 前 的 竖 向 位 移 ;2 代表 改进 后 的 竖 回 位 移 。 
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(a) 沿 图 11 上 部 所 示 (b) 沿 图 11 下 部 所 示 


图 15 沿 纵 筋 轴 疝 的 各 截面 最 大 正 应 力 对 比 


Fig. 1$ Comparison of maximum normal stress along the longitudinal reinforcement 


投稿 网 站 :http://cjam. xjtu. edu. cn 微 信 公众 号 :应 用 力学 学 报 


402 学 报 第 41 卷 
50 本 二 本 、 
i ee 改进 前 一 一 改进 后 
和 2 
= = 
RR RR 
武 武 
到 到 
监 鉴 
和 和 
0 5 10 1 20 0 5 10 15 20 
时 间 /s 时 间 /s 
(a) 2 所 桥墩 (b) 3 与 桥墩 
图 16 墩 确 筋 力 时 程 对 比 
Fig. 16 Comparison of pier-bottom shear time history 
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图 17 党 桥墩 高 度 的 最 大 剪 力 对 比 
Fig.17 Comparison of maximum shear along the piers 
= 日 上 
[1] 周 田 宇 , 梁 富源 , 马 完 君 ,等 . 川 演 区 域 高 频 地 震动 地 形 效 应 初 


窗 水 峡谷 -丘陵 地 貌 地 震动 的 合理 输入 ,是 开展 
坐落 于 此 类 地 形 长 大 工程 抗震 研究 的 基本 前 提 。 本 
研究 针对 此 研究 空白 与 现实 需求 , 围 比 相 干 函 数 模 
型 的 建立 .模型 自身 特征 和 对 地 震动 影响 ,以 及 对 结 
构 地 震 反 应 的 影响 等 三 方面 ,依次 开展 研究 。 具 体 
得 出 如 下 结论 。 

1 ) 提 出 并 推导 得 到 了 和 窗 水 峡谷 -丘陵 复合 地 貌 
相干 兄 数 的 表达 式 , 同 时 引入 了 撕 陵 土 体 的 滤波 效 
应 和 竖 向 距离 对 相干 性 的 有 影响。 与 改进 前 模型 相 
比 , 相 干 性 偏 低 , 原 因 在 于 竖 向 距离 增加 了 传播 路 径 
的 差异 ,进而 幅 值 和 相位 角 发 生变 异 概 率 增 加 。 

2 ) 编制 程序 模拟 了 基于 本 人 研究 模型 的 峡谷 -于 
发 复合 地 貌 地 下 多 点 地 震动 ,验证 了 其 合理 性 ,对 比 
结果 发 现 : 改 进 后 不 同 目标 点 间 地 震动 标准 差 间 的 
差 值 增 大 ,表明 基于 本 研究 相干 函数 的 各 点 地 震动 
的 空间 变异 性 明显 增强 。 

3 ) 为 进一步 表明 本 研究 意义 ,对 一 刚 构 桥 地 震 
反应 进行 分 析 对 比 , 结 果 表 明 : 基 于 本 研究 模型 的 地 
震动 输入 下 的 结构 地 震 差 动 反 应 也 明显 增 大 ,说 明 
若 忽 略 复合 地 形 地 貌 对 相干 函数 影响 ,会 错 佑 甚至 
低估 结构 地 震 反 应 ,从 而 导致 安全 隐患 。 
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